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ABSTRACT 

The bright-rimmed cloud SFO 22 was observed with the 45 m telescope of Nobeyama Radio Ob- 
servatory in the 12 CO (J = 1-0), 13 CO (J = 1-0), and C 18 (J = 1-0) lines, where well-developed 
head-tail structure and small line widths were found. Such features were predicted by radiation-driven 
implosion models, suggesting that SFO 22 may be in a quasi-stationary equilibrium state. We com- 
pare the observed properties with those from numerical models of a photo-evaporating cloud, which 
include effects of magnetic pressure and heating due to strong far-ultraviolet (FUV) radiation from 
an exciting star. The magnetic pressure may play a more important role in the density structures of 
bright-rimmed clouds, than the thermal pressure that is enhanced by the FUV radiation. The FUV 
radiation can heat the cloud surface to near 30 K, however, its effect is not enough to reproduce 
the observed density structure of SFO 22. An initial magnetic field of 5 fj,G in our numerical mod- 
els produces the best agreement with the observations, and its direction can affect the structures of 
bright-rimmed clouds. 

Subject headings: ISM: clouds — ISM: magnetic fields - ISM: individual (SFO 22) - Methods: numer- 
ical - Methods: observational 
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1. INTRODUCTION 

Bright-rimmed clouds (BRCs) are cometary molecu- 
lar clouds found at the edge of H II regions. These 
clouds have bright rims on the side facing the excit- 
ing star and extended tails on the other side. Since 
their head-tail morphologies suggest that BRCs are in- 
teracting with the radiation or stellar wind from an ex- 
citing star, BRCs are considered as potential sites of 
triggered star formation. The gradients of age spread 
in young stars along the axes of BRCs indicate that 
their formatio n may have been sequentially triggered by 
shock waves (iSugitani et al.lll995t IGetman et al.l 120071 : 
llkeda et alj|2008t IGetman et al.ll2009fi ~ 

Radiation-driven implosion models are often consid- 
ered for the formation and evolution of BRCs and their 
triggered origins. Strong UV radiation from nearby mas- 
sive stars can photoionize and photoevaporate the sur- 
faces of surrounding molecular clouds, and whose ef- 
fects have been studied by various groups. IBertoldil 
(1989) developed an approximate analytical solution 
for the evolution of molec ular cloud compres s ed by 
radiation-driven implosion. iLefioch k, Lazarefll (1994) 
investigated a radiation-driven implosion model us- 
ing hydrodynamic simulations. Recent hydrody- 
namic simulations of radiation-driven implosion in- 
clude effects of physics such as self-gravity of the gas 
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fuse radiation field (jHaworth fc Harries 2012 ), and tur - 
bulence in molecular cl ouds (jGritschneder et al.l 120091) . 
IMotovama et al.l ((2007) demonstrated that radiation- 
driven implosion can enhance accretion rates enough to 
account fo r the high luminosit ies of YSOs observed in 
the BRCs (jSugitani et al.lll989l ). Typical shock speed of 
a few km s _1 in this study is consistent with that esti- 
mated from observat ions of age gradient s of yo ung stars 
in and aroud BRCs (IGetman et al.ll2007l 120091 ). 

Radiation-MHD studies suggest the possibility of al- 
tering evolution of p hotoionized cl o ud by the presence 
of magnetic field. IHennev et al.l (|2009l ) carried out 
the first three-dimensional radiation-MHD simulations 
of photoionization of a magnetized molecular globule 
under ultraviolet radiation. They reported the photo- 
evaporating globule will evolve into a more flatten shape 
compared with the non-magnetic case when the cloud ini- 
tially has a strong magnetic field (that is 100 times the 
thermal pressu r e) pe rpendicular to UV radiation field. 
IMackev fc Liml (|201 ID also showed that strong magnetic 
field has significant influence on the dynamics of the pho- 
toionization process. Measuring magnetic field strengths 
in BRCs observationally has been difficult, and the ef- 
fects of magnetic field can only be inferred indirectly by 
comparing the observed density structures and kinemat- 
ics with those obtained by theoretical models. 

Strong far-ultraviolet (FUV) radiation from an ex- 
citing star may also influence the evolution of a pho- 
toionized cloud. As shown in many studies of pho- 
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ton dominated regions (e.g. iTielens k HoUenbachl 119851: 
Holl enbach et al.lll991l) . the FUV radiation from a mas- 
sive star can heat molecular clouds through photoelectric 
heating and photodissociation of important coolants such 
as carbon monoxide. Temperature of the cloud heated 
by FUV radiation ranges from a few tens up to a few 
hundreds of Kelvin, strongly dependent on the intensity 
of the FUV radiation and the density of cloud. High 
thermal pressure enhanced by FUV radiation may affect 
the evolution of BRCs. However, there have been lim- 
ited theoretical works that i nclude heating due to FUV 
radiation (jMiao et al.ll2006t iHennev et al.l l2009h . These 
effects should be included in theoretical models for reli- 
able comparisons with observations. 

The purpose of this study is to investigate how mag- 
netic field and FUV radiation affect the evolution of 
BRCs through the comparisons between observations 
and numerical models. An evolutionary scenario of BRCs 
due to radiation-driven implosion has an initial implosion 
phase followed by a quasi-stationary equilibrium phase. 
In this paper, we focus on the quasi-stationary equilib- 
rium phase, and implosion phase will be investigated in a 
subsequent paper. The BRC SFO 22 was observed with 
the 45 m telescope of Nobeyama Radio Observatory, and 
the results were compared with numerical models of pho- 
toionized clouds with effects of magnetic field and FUV 
radiation. The layout of this paper is as follows. In 
section the details of observations are described. In 
section [3j a description of our numerical models is given. 
Section H] and section [5] give the results and discussions. 
In section [51 we summarize our main conclusions. 

2. OBSERVATIONS AND ANALYSIS 
2.1. Observed Bright-rimmed Cloud 

BRC SFO 22 is located at the eastern edge of H II 
region s281, and w as selected from BRC catalog of 
ISugitani et al.l (|199lD . Figure Q] (a) shows the entire im- 
age of the region s281, andQ](b) gives the close-up view 
of SFO 22 enlarged from[T](a). This H II region is ionized 
by 6' 1 Ori, marked by a cross in Figure[T](a). The spectral 
type of the primary exciting star in 1 Ori is 07V, and has 
a pro jected distance of 6.5 pc to SFO 22 (| Morgan et al.l 
200 4. S281 (lBlaauw|[r9 64) is 460 pc away from us. 

ISugitani et al.l (|1991| ) classified the BRCs into three 
types in order of increasing degree of rim curvature: type 
A, B, and C. Based on the radiation-driven implosion 
model, the shapes of cloud rims reflect the evolution- 
ary stages. The type A BRCs are still undergoing com- 
pression by the Shockwaves generated by ionization, and 
type B and type C BRCs are approaching or have al- 
ready reached the phase of quasi-stationary equilibria. 
SFO 22 has a well-developed head-tail structure along 
the line to the exciting star and it is classified as type B. 
Although the head part of the cloud contains the IRAS 
point source 05359-0515, ammonia rotational inversion 
lines were not detec ted toward this IRAS point source 
([Morgan et all feOlO). Since ammonia lines trace dense 
gas associated with protostellar cores, the non-detection 
of ammonia lines has been interpreted as no star forming 
activities in SFO 22. 

2.2. Observations 
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Fig. 1. — (a) DSS red image of the H II region s281. The cross 
labels the position of the exciting star, (b) DSS red image of bright- 
rimmed cloud SFO 22. The arrow indicates the direction to the 
exciting star. 

(A color version of this figure is available in the online journal.) 

Observations were carried out using the 45 m telescope 
of Nobeyama Radio Observatory in 2005 January and 
April. We observed 12 CO (J = 1-0) at 115.271203 GHz, 
13 CO (J = 1-0) at 110.201370 GHz, and C 18 (J = 1- 
0) at 109.782182 GHz. The half-power beam width of 
the telescope and main-beam efficiency at 115 GHz were 
15"and tjmb = 0.4, respectively. We used the 5x5 beam 
focal plane array receiver "BEARS" whose beam sepa- 
ration is 41.1". As receiver backends we used a 1024 
channel digital autocorrelator with a 31.25 kHz frequency 
resolution. The corresponding velocity resolution is 81.5 
ms" 1 at 115 GHz. The typical system noise tempera- 
ture was 300-450 K depending on the atmospheric con- 
ditions. The intensity scale of the spectra was calibrated 
by the chopper wheel method. The corrected antenna 
temperature T* A is converted into main beam brightness 
temperature using the relation of Tg = T^/r] m b. 

We observed SFO 22 with a grid spacing of 20.55", 
which is half of the beam separation of BEARS, in the 
line of 12 CO (J = 1-0). The mapped area was 390" x 
390". Dense regions of clouds were observed with finer 
grid spacing of 10.3"in the lines of 13 CO (J = 1-0) and 
C 18 (J = TO). The mapped area for these lines were 
195" x 195" and 154" x 154", respectively. The pointing 
accuracy of the antenna was checked and corrected every 
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1.5-2 hr using SiO maser emission from Ori KL, and its 
typical error was less than 5". All observations were 
carried out by position switching mode. The data were 
reduced by using the software package NewStar provided 
by Nobeyama Radio Observatory. 

2.3. Observational Results and Analysis 

Figure [5] shows the velocity-integrated intensity maps 
of SFO 22. The reference center of the map is the peak 
position of C 18 (J = 1-0) emission at RA(2000) = 5 h 
38 m 21.6 s , Dec(2000) = -5° 13' 37.8". Emission from 
12 CO (J = 1-0) coincides with the optical images. The 
cometary morphology is clearly shown. On the contrary, 
the C 18 (J = 1-0) emission is very weak and detected 
only at a few points. 

Figure [3] shows the observed line spectra toward the 
peak position of C 18 (J = 1-0) and offset positions along 
declination. Assuming that 12 CO (J = 1-0) emission is 
optically thick at the peak of C 18 (J = 1-0) emission, 
we can calculate the excitation temperature as 



5.53 



ln[l + 5.53/(T B ( 12 CO) + 0.819)] ' 



(1) 



where Tb( 12 CO) is the brightness temperature of 12 CO 
(J = 1-0) at the peak of C 18 (J = 1-0). The excita- 
tion temperature of SFO 22 is found to be 27.1±1.8 K. 
Fig. S] shows the position- velocity diagrams of 12 CO (J 
= 1-0), 13 CO (J = 1-0), and C 18 (J = 1-0) along dec- 
lination through the peak position of C 18 (J = 1-0). 
The line widths of 12 CO (J = 1-0), 13 CO (J = 1-0), and 
C 18 (J = 1-0) are roughly - 2.0 km s" 1 , - 1.5 km 
s _1 , and ~ 1.0 km s , respectively. Line widths of SFO 
22 are relatively narrow compared to other BRCs. In 
the millime ter and sub-mil l imete r molecular line survey 
of BRCs bv lDe Vries etHl jUSI), many BRCs have CO 
line widths of > 5 km s _1 . These large line widths of 
BRCs may be attributed to large velocity dispersion due 
to turbulence and outflow activities. The observed nar- 
row line widths of SFO 22 could suggest that the influ- 
ence from the photoionization-induced shocks may have 
already disappeared, and SFO 22 has now reached the 
phase of quasi-stationary equilibrium predicted in the ra- 
diation implosion model. 

Figure [5] shows column density profiles of SFO 22 along 
declination, which is almost parallel to the direction to 
the exciting star. We derive column density distribution 
of the cloud using data from this observation. Assuming 
that the excitation temperatures of the 13 CO and C 18 
lines are the same as the 12 CO line, we can calculate 
optical depth of 13 CO and C 18 using 
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T s ( 13 CO) 



5.29(l/(exp(5.29/T ex ) - 1)) - 0.164) 



(2) 



Tb(C 18 0) 



5.27(l/(exp(5.27/T e:c ) - 1)) - 0.166) 



(3) 

respectively. The column densities of 13 CO and C 18 
molecules can be derived as 



A( 13 CO) = 2.42 x 10 14 r Tl3Awi3T - 



exp[-5.29/T e J 



(4) 
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Fig. 2.— Integrated intensity maps of SFO 22. The 12 CO (J = 
1-0) map (top) has a lowest contour of 1.39 K km s _1 (3 a) and 
contour intervals of 1.39 K km s" 1 (3 a). The 13 CO (J = 1-0) 
map (middle) has a lowest contour of 0.315 K km s — 1 (3 tr) and 
contour intervals of 0.525 K km s" 1 (5 cr). The C ls O (J = 1-0) 
map (bottom) has a lowest contour of 0.131 K km s (3 a) and 
contour intervals of 0.0435 K km s _1 (1 a). 
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Fig. 3.— Some selected 12 CO (J = 1-0), 13 CO (J = 1-0), and C ls O (J = 1-0) spectra observed toward the SFO 22 along declination 
through the peak of C ls O (J = 1-0). 
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Fig. 4.— Position-velocity diagrams of 12 CO (J = 1-0) (left), 13 CO (J = 1-0) (center), and C ls O (J = 1-0) (right) along declination 
through the peak positions of C 18 (J = 1-0) emission. The diagram for 12 CO (J = 1-0) has lowest contour of 1.53 K (3 <r) and contour 
interval of 1.53 K (3 a). The diagram for 13 CO (J = 1-0) has lowest contour of 0.276 K (3 <r) and contour interval of 0.92 K (10 a). The 
diagram for C 18 (J = 1-0) has lowest contour of 0.125 K (3 a) and contour interval of 0.125 K km (3 a). The vertical dotted lines indicate 
the systemic velocity of 10.4 km/s. 
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iV(C iB 0) = 2.42 x 10 1 



ri 8 Awi 8 r e . 



1 - exp[-5.27/T e: 



(5) 



where Avi3 and Auig are line widths of the 13 CO 
(J = 1-0) and C 18 (J = 1-0) emission, respectively. 
The column density of 13 CO is converted to the col- 
umn density of H2 by assuming an abundance ratio of 
N(H 2 )/iV( 13 CO) = 5.0 x 10 5 (jDickmanl Il97cl . Un- 



der the environment where molecular clouds are illu- 
minated by strong UV radiation, abundance ratio of 
7V(H 2 )/iV(C 18 0) are thought to be affected by sele ctive 
destr uction by UV radia tion (|Glassgold et al.lll985[ ). We 
follow [Ni wa et all (j2009) to obtain the column density of 
H2 from the column density of C 18 0. The column densi- 
ties of 13 CO and C 18 can be formulated by least-square 
fitting as 



iV(C 18 0) = 3.91 x 10~ 2 x 7V( 13 CO). 



(0) 
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Fig. 5. — Column density profiles of SFO 22 along the declina- 
tion. The open circles (colored blue in the online version) and filled 
circles (colored red in the online version) label the column densi- 
ties derived from 13 CO (J = 1-0) emission and C ls O (J = 1-0) 
emission, respectively. 

(A color version of this figure is available in the online journal.) 



An abundance ratio of N(R 2 ) /N(C ls O) = 1.28 x 10 7 can 
be derived. This value is slightly larger than the standard 
value of 6.0 x 10 6 for molecular clou ds not associated with 
H II regions (|Frerking et al.|[l982f) . Figure [5] shows that 
the column density profiles obtained here are nearly flat 
with the column density of ~ 10 22 cm" 2 . Adopting a 
distance of 460 pc to SF022, the total masses traced by 
13 CO emission and C 18 emission are M i3 = 12.0 M Q 
and Mig = 5.3 M , respectively. 

3. MODEL DESCRIPTION 

Our numerical model is based on the analytic model of 
IBertoldi fc McKeel (|1990D for a photo-evaporating cloud 
in the quasi-stationary equilibrium state with a poly- 
tropic equation of state. The polytropic gas is a good 
approximation for the cases where either gas pressure or 
magnetic pressure is dominant. In an actual bright- rimed 
cloud, these two quantities may have been comparable 
before the compression by radiation-driven implosion al- 
ter the state, and FUV radiation will affect the thermal 
structure through photoelectric heating and photodisso- 
ciation of molecular coolants. Here we adopt a more re- 
alistic model with the inclusion of thermal pressure and 
magnetic pressure explicitly and the heating due to FUV 
radiation. 

3.1. Density Structure 
Total pressure of the gas is expressed as 



Pint. — Pi 



+ Pr 

7 8ir ' 



C 2 S P 



(7) 



where, P t h, Pmagi c s, P, 7, and B are the thermal pres- 
sure, the magnetic pressure, the sound speed, the density 
of the cloud, the ratio of specific heats, and the mag- 
netic field strength, respectively. A 7 = 5/3 is adopted, 
which is appropriate for the molecular clouds because 
the molecular hydrogen behaves like a monoatomic gas 
at temperature < 100 K. We assume ideal gas, so that 
the sound speed c s is related with the gas temperature 
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Fig. 6. — Schematic figures of compression of magnetized cloud 
due to radiation-driven implosion for two extreme cases. UV ra- 
diation propagates downward, and the cloud is compressed only 
along this direction in both cases, (left) The magnetic field is par- 
allel to the direction of UV radiation, (right) The magnetic field is 
perpendicular to the direction of UV radiation. 



T as 



(8) 



where ks and fi are the Boltzmann constant and mean 
mass per nucleus in unit of the hydrogen mass mu = 
1.67 x 10~ 24 g, respectively, and fi = 1.15. For simplicity, 
the effects of magnetic field is approximated through 



(9) 



B = B n \ — 

Pa 



where Bq and po are magnetic field strength and density 
of the cloud before undergoing compression by radiation- 
driven implosion. The exponent a, which ranges from 
to 1, is a parameter representing how much magnetic field 
is trapped in the gas during the compression. Figure [6] 
shows schematic drawings of two extreme cases by as- 
suming that compression of a cloud perpendicular to the 
direction of radiation is small. If the cloud is compressed 
along the magnetic field, the magnetic field strength will 
hardly change during compression, i.e. a ~ 0. On the 
other hand, if the cloud is compressed perpendicular to 
the magnetic field, the strength of the field increases as 
B (x p during compression, i.e. a ~ 1. The value of a 
depends on the initial configuration of magnetic field and 
the shape of cloud. We leave a an open parameter as it 
is hard to determine an accurate value without launch- 
ing MHD simulations. We also neglect the diffusion of 
magnetic field due to the longer timescale of ambipo- 
lar diffusion compared to the dynamical timescale of the 
radiation-driven implosion. 

We calculate structure of the cloud in axisymmetry 
along the line to the exciting star. We also assumed 
that distance from the cloud to the exciting star is larger 
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Fig. 7. — The coordinate system used in our numerical models. 
The exciting star is assumed far above. Ionizing UV radiation is 
parallel to the symmetric axis. The flow of evaporation streams 
along the surface normal. The angle between surface normal and 
UV radiation is denoted by 9. 



than size of the cloud, so that UV radiation field can 
be treated as planar. Figure [7] illustrates the coordi- 
nates system we use in this paper. UV radiation prop- 
agates downward parallel to the z-axis, and the ionized 
gas e vaporates off the c l oud s urface with angle 9 to z- 
axis. lBertoldi~fc McKcc (1990) showed that the position 
of the cloud surface can be approximated as 



R c a 2 lncos(ar/i? c ), 



(10) 



where R c and a are the curvature r adius at z = and 
the clo ud width parameter defined in IBertoldi fc McKeel 
(1990), respectively. The method used to determine a 
and R c is described in section [5751 

In this study, we assume that all part of the cloud is 
subject to the same acceleration through the rocket effect 
of the photo-evaporation flow, as the cloud is in quasi- 
stationary equilibrium state. The equation of hydrostatic 
equilibrium for the cloud is 



dP t , 



dz 



gp, 



(ii) 



where g is the acceleration of the cloud. With the equa- 
tion (0 and (J9)), this equation can be rewritten as 



dp _ gp - 2c s p/~f(dc s /dz) 

dz ~ c 2 / 7 + 2af? 2 0/ p ) 2a -78^o ' 



(12) 



The density distribution inside the cloud is determined 
by solving this equation with appropriate boundary con- 
ditions. 

The boundary conditions at cloud surface are deter- 
mined using jump conditions for D-critical ionization 
front. Therefore, pressure at cloud surface is given by 



Psur(r) = 2pmnCiF uv (r) cos 9, 



(13) 



where Ci = lOkms -1 and Fjjv( r ) are the sound speed of 
ionized gas and ionization photon flux reaching ionization 
front, respectively, and the width of the ionization front 
is negligible. Since some part of the incident ionizing 
UV photons is consumed by recombined hydrogens in the 
photoevaporation flow, ionizing UV photon flux arriving 
at the cloud surface is written as 



a B n e (r, z)n p (r, z) dz, 



(14) 



where Fi, olb = 2.7 x 10~ 13 cm s , n e , and n p are the 
incident ionizing photon flux, the hydrogen electronic re- 
combination coefficient into the excited state, the elec- 
tron number density, and the proton number de nsity, re- 
spectively. Following IBertoldi fc McKeel (|1990l ). we in- 
troduce a dimensionless parameter lo that represents the 
effective fractional thickness of the recombination layer. 
Equation (|14l) can be rewritten as 



Fuv(r) = Fi — Lu{r)a B n 2 n (r = 0)R C , 



(15) 



where nn is the hydrogen number density just behind 
the ionization front. If we assume the ionized gas streams 
away from the ionization front with the velocity of sound 
speed, and the stream line is normal to the cloud surface, 
we can derive the form, 



w(r)~3w(0) 



f 


\r'P to t(r')] 


Jo 


I rPtot(0) \ 




R(r>) 


sin( 


V + r - r'] 



1 2 



dr' 



(16) 



where R(r') is the curvature radius of the cloud surface 
at r = r'. We adopt u> for a spherical cloud with radius 
R c at the cloud tip: 



w(0) 



g(g-l) 



(17) 



where q and ip are the ratio of incident ionizing photon 
flux Fi to ionizing photon flux reaching ionization front 
Fuv and the photoevaporation parameter, respectively. 
The photoevaporation parameter is defined as 



asFiRc 



(18) 



The value of ip in our numerical models ra nges from 74 to 
104 depending on the model parameters. iSpitzerl (|1978l ) 
derived an analytic estimation of q as 



q 



1 + (1 + 1.5^ 1/4 ) : 



(19) 



The density and pressure distributions of the cloud are 
characterized by scale height defined as h c = Ptot{z — 
0)/gp(z = 0). If i/j > 1, which is the case for SFO 22, h c 
is related to R c as R c — Q.5h c . 

3.2. Thermal and Chemical Model 

We solve the reaction networks for the species of H 2 , 
CO, C + , O, and the electron e. We adopt the simplified 
reaction model as described in iNelson fc LangeH (1997) 
to determine the abundance of CO molecules. In this 
reaction model, C + is directly converted to CO without 
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accounting explicitly for the intermediate reactions. The 
formation rate of CO molecules is expressed as 



Rco = 5 x 10" ib n(C+)n(H 2 )/3, 



(20) 



where n(C + ) and n(H2) are the number densities of C^ 
and H2, respectively. The coefficient f3 is defined as 



5 x l(T 10 X(O) 



5 x 10- 10 X(O) + £>ch x /«(H 2 ) ' 



(21) 



where A(O) and -Dch x are the fractional abundance of 
oxygen and the total photodissociation rate of both CH 
and CH2. This total photodissociation rate is written as 



D 



CH X 



= 5 x 10- iu G exp(-r t/ y) 



(22) 



where Go and t\jv are the intensity of the incident FUV 
radia tion in terms of the Habing interstellar radiation 
field ()Habing||1968ft and the optical depth in UV range, 
respectively. The optical depth tjjv is related to the 
visual extinction A v by tjjv = 2.5A V . We adopt the 
conversion factor between the visual extinction and the 
total hydrogen column density as Xa v — Ay/Nn = 6.3 x 
10 -22 mag cm 2 in this paper. The photodissociation rate 
of CO by the FUV radiation is written as 



.Deo = 10- 10 G exp(-T C/y )n(CO), 



(23) 



where n(CO) is the number density of CO. The abun- 
dance of CO is calculated by solving the equation of for- 
mation and dissociation balance 



Rco - D co - 0. 



(24) 



Abundances of other species are calculated as follows. 
The cloud is assumed to be composed of molecular hy- 
drogen, so that the number density of H2 is written as 



n(H 2 ) = 0.5n, 



(25) 



where n is the number density of hydrogen nuclei and 
related with mass density as 



P 



(26) 



We assume that carbon exists in ionized form in the cloud 
as carbon is easily photoionized by FUV radiation ow- 
ing to its lower ionization energy (11.2 eV) compared to 
that of hydrogen. Although this assumption may cause 
overestimation of cooling rate through C + in the deep 
inner region of cloud where FUV is strongly attenuated 
and rotational line emission from CO is a more impor- 
tant cooling process. Our assumption docs not affect the 
numerical results, however. The number density of C + 
is calculated by 



n(C+) = X{C tot )n - n(CO). 



(27) 



where the elemental abundance of carbon is taken to be 
V(Ctot) = 10~ 4 . We assume that oxygen exists in atomic 
form in the cloud, because its ionization energy (13.6 eV) 
is similar to that of hydrogen. The number density of 
oxygen is calculated by 



n(0) = X(O to t)n-n(CO). 



(28) 



where the elemental abundance of oxygen is taken to 
be A(Otot) = 2.0 x 10~ 4 . Constant electron fraction 



n e /n — 10~ 7 is also assumed to calculate number density 
of electron. 

For the heating processes in the cloud, we consider 
photoelectric heating and cosmic ray heating. T he pho- 
toelectric heating rate is (jBakes fc Tielenslll994f) 

r pe = 10 _24 eG exp(-r; 7 y)rt ergs cm~ 3 s -1 , (29) 

where e is the photoelectric heating efficiency 

4.87 x 10~ 2 



1 + 4 x 10- 3 (G exp(-r c/ y)T 1 /2/n e ) - 73 

3.65 x 10- 2 (r/10 4 ) - 7 
+ 1 + 2 x 10- 4 (G exp(-T l /y)T 1 /2/ ne ) ' 



(30) 



Cosmic ray heating become the dominant heating pro- 
cess in the inner region where FUV radiation does not 
penetrate. Cosmic ray heating rate is given by 



( p (H 2 )AC; cr .n(H2) ergs cm 3 s~ 



(31) 



where C P (H 2 ) is the primary cosmic ray ionization rate of 
H 2 and AQ cr is the energy deposited as heat as a result 
of this ionization. We adopt values of Cpi^) — 7.0 x 
10~ 17 s" 1 (Ivan Dishoeck fe Black| [T986h and AQ cr = 20 
eV (IGoldsmith fc Langerlll978ft~ ~^ 

For the cooling processes in the cloud, we consider ra- 
diation from the CO rotational transitions, collisionally 
excited line emission from C + and O, and collisional heat 
transfer between gas and dust. The cooling rate due to 
CO, Aco is t aken from the tabul ated cooling function 
computed by INeufeld et al.l (fT995ft for T < 100 K and 
INeufeld fc Kaufman! (11993ft for T > 100 K. The cooling 
rates due to collisional excitat ion of C + and O are taken 
from iNelson fc Langerl (|1997ft . The cooling rate due to 
C+ is 

A c + = [1.1 x 10- 23 n(C+)ncxp(-92/T)] 

/{I + HH 2 )K„ t )[l + 2exp(-92/T)]}, (32) 

where the critical number density is taken to be n cr it = 
3 x 10 3 cm~ 3 . The cooling rate due to O is 



A = 5.0 x 10~ 27 n(O)n[24exp(-228/T)] 



+7exp(-326/T)]T 



1/2 



(33) 



The cooling rate by the dust grain is 
(|Hollenbach fc McKed [r989) 



A dust = [1.2 x 10 



" 31 n 2 



/ T 



1/2 



1/2 



ziooAy 
Vioook; ^ a mm J 

x [1 - 0.8 exp(-75/T)] (T - T dust ),(34) 

where <z m j n and Tdust are the minimum radius of grains 
and the dust temperature, respectively. We used a m i n — 
100 A to calculate the cooling rate by dust grains. The 
dust temperature is cal culated following the method by 
IHollenbach etaLl (fl99l . 

3.3. An Iterative Procedure for Obtaining Numerical 

Solution 

An iterative procedure is used to achieve the density 
distribution of the cloud. We use a uniform grid of 
200 (radial) x 1000 (axial) cells. The grid spacing is 
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6.5 x 10 4 pc. For presentation, we symmetrize the re- 
sults obtained on a grid of 400 x 1000 cells. In order 
to fit the width of the cloud with that of SFO 22, the 
cloud width parameter a is calculated by substituting 
r = 0.13 pc and z — —0.6 pc into equation ([TUl) . The 
procedure is as follows, (i) Equation (|12j) is numerically 
solved by Runge-Kutta method using current tempera- 
ture and sound speed, (ii) Optical depth to the cloud 
surface is calculated, (iii) Chemical reaction network 
is solved to determine abundances of included species, 
(iv) The temperature and sound speed are updated us- 
ing new chemical abundances. The thermal equilibrium 
is assumed to determine the temperature of cloud: 



pe 



Acq — A c + — Aq — Adust — 0. 



(35) 



These processes are repeated until the converged solu- 
tion is achieved. We have confirmed that relative errors 
between the two consecutive steps are smaller than 10~ 2 
for all models in this paper. 

Once the curvature radius of the cloud R c is given, we 
can determine the position of cloud surface from equation 
(fT0|) and calculate structure of the cloud using iterative 
procedure described above. Since we assume that the 
cloud is in an equilibrium state, we adopt R c which min- 
imizes the pressure gradient along r-direction. We change 
R c from 0.02 pc to 0.11 pc in increments of 0.001 pc, and 
calculate structure of the cloud for each R c . Then, we 
integrate deviation of total pressure Ptot (t, z) from the 
average total pressure over all cells at the same z coordi- 
nates P ave (z), 



j \(P to t(r, z) - P ave {z)) /P ave (z)\dV 
JdV 



(36) 



We adopt the structure of cloud for which this value is 
minimized. 

4. RESULTS 

We present results of our numerical models in this 
section. Model parameters of the run are summarized 
in Table [T] To compare with our observational re- 
sults, intensities of the incident ionizing UV and FUV 
radiation were determined to correspond to those ex- 
pected in the region where SFO 22 is located. Since 
the spectral type of the primary exciting star of SFO 
22 is 07V, we adopted the values of log Suv = 48.76 
s- 1 and log Sfuv = 48.76 s" 1 as the UV and FUV 
photo n luminosities (|Vacca et alJll996l iDiaz-Miller et all 
1998). These luminosities and the distance from the ex- 
citing star to SFO 22 of 6.5 pc give the incident ion- 
izing UV flux of Fi — 1.192 x 10 9 cm~ 2 s _1 and inci- 
dent FUV flux of Go = 9 4 in terms of the Habing field 
F H = 1.21 x 10 7 cm" 2 s" 1 (jBertoldi fc Drainelll996l i. We 
used these values in all models except model A. The 
number density of the cloud before undergoing compres- 
sion by radiation-driven implosion was assumed to be 
no = po/ limn = 10 3 cm -3 in all models. Some physi- 
cal values obtained from our results are summarized in 
Table [2] In Figure [H column density profiles along the 
z-axis for all models are plotted with column densities of 
SFO 22 derived from our observations. 

4.1. Effects of Strong FUV Radiation and Magnetic 

Field 



TABLE 1 
Model parameters 



Model 


a 


B [flG] 


Go 


Fi [cm ' l s i ] 


A 


- 





1.0 


1.192 x 10 u 


B 


- 





94.0 


1.192 x 10 9 


CI 


0.25 


25 


94.0 


1.192 x 10 9 


C2 


0.25 


45 


94.0 


1.192 x 10 9 


C3 


0.25 


60 


94.0 


1.192 x 10 9 


C4 


0.25 


80 


94.0 


1.192 x 10 9 


Dl 


0.50 


7 


94.0 


1.192 x 10 9 


D2 


0.50 


15 


94.0 


1.192 x 10 9 


D3 


0.50 


25 


94.0 


1.192 x 10 9 


Dl 


0.50 


15 


94.0 


1.192 x 10 9 


El 


0.75 


2 


94.0 


1.192 x 10 9 


E2 


0.75 


5 


94.0 


1.192 x 10 9 


E3 


0.75 


10 


94.0 


1.192 x 10 9 


E4 


0.75 


25 


94.0 


1.192 x 10 9 
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Fig. 8. — Comparisons of the column density profiles along the z- 
axis with the column densities derived from observations. The top, 
middle, and bottom panels show the results with a = 0.25, 0.50, 
and 0.75, respectively. The dashed line (colored green in the online 
version) represents the column density profile of model A. Column 
densities of SFO 22 are plotted with same symbols as Figure [5] 
(A color version of this figure is available in the online journal.) 
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In this subsection, we describe results of following three 
typical models, A, B, and E2 to show how strong FUV 
radiation and magnetic field affect structures of clouds. 
Figure [9] shows the densities, the column densities, and 
the temperatures of these models. In model A, as a ref- 
erence, we calculated the structure of the cloud without 
magnetic field, assuming strength of average interstellar 
FUV radiation field of Go = 1. In model B, to see the 
effects of strong FUV radiation, we calculated the struc- 
ture of the cloud assuming strength of FUV radiation 
expected at the region where SFO 22 is located, but the 
effects of magnetic field was not included. In model E2, 
we included not only the effects of strong FUV radiation 
but also the effects of magnetic field. The initial mag- 
netic field strength Bq and a were set to be 5 pG and 
0.75, respectively. 

Figure [3] shows that the column density distribution of 
model A is much different from that of observed cloud. 
The peak column density of 1.26 x 10 23 cm -2 is one order 
of magnitude higher than that of observed cloud (~ 2 x 
10 22 cm -2 ). In addition, slope of column density profile 
is steeper than observations. Figure [§] shows that the 
cloud has nearly constant temperature of 20 K in model 
A. Since FUV radiation is not significant in this model, 
effective cooling keeps the dense region relatively cold. 

Comparisons of model B to model A show that FUV 
radiation has little influence on the density and the ther- 
mal structures of the cloud. Although the cloud surface 
is heated near 30 K in model B, this warm surface layer is 
very thin. High density at the head region prevents CO 
molecules, which are main coolant of molecular gas, from 
photodissociation except thin surface layer. Isothermal 
gas is good approximation at inner region of the cloud. 
Figure [5] shows that the column density of model B is 
much higher than observations as well as model A. Al- 
though FUV radiation reduces the density than model A 
by factor of a few, it is not enough to reproduce observed 
low column densities. The shape of the cloud slightly dif- 
fers from model A. Strong FUV radiation enhances ther- 
mal pressure at head region and makes the curvature 
radius at cloud tip larger than that of model A. Moving 
to tail side, the differences of the density and the column 
density from those of model A become larger. As a re- 
sult, the slope of the column density profile of model B 
is steeper than that of model A. 

Comparisons of model E2 to model A and model B 
show that magnetic field reduces the density of the cloud. 
Figure [9] shows that the density and column density of 
model E2 at head region are one order of magnitude lower 
than those of model A and B. As shown in l4.2[ magnetic 
pressure is dominant at z > —0.4 pc in model E2. Ad- 
ditional support due to magnetic pressure makes density 
of the cloud lower than models without magnetic field. 
The warm surface layer of model E2 is thicker than that 
of model B, because lower density at head region allows 
FUV radiation penetrating deeper inside the cloud. Fig- 
ure|8]shows that the slope of the column density profile of 
model E2 is flatter than those of model A and B. The col- 
umn density profile of model E2 shows better agreement 
with that of the observed cloud than other two models. 

4.2. Dependence on the Density Dependence of 
Magnetic Field Strength 



In this subsection, we present comparison of following 
three models, C2, D2, and E2 to show how the value 
of exponent a in equation ([9]) affects the results. The 
values of a and Bo in model C2, D2, and E2 were set to 
be 0.25 and 45 fiG, 0.50 and 15 fiG, and 0.75 and 5 /xG, 
respectively. Comparison of these models reveals that 
the value of a influences structure of cloud. 

As can be seen from Figure HI although these three 
models have similar maximum values of column den- 
sity of ~ 2 x 10 22 cm~ 2 , the column density profiles of 
these models are qualitatively different. The model with 
smaller value of a has steeper column density profile at 
the head region (z > -0.2 pc) and has flatter column den- 
sity profile at tail region (z < -0.3 pc). In model C2, the 
slope of the column density profile becomes flatter with 
moving to tail side. Opposite trend is observed in model 
E2. The slope of the column density profile becomes 
steeper with moving to tail side. Contrary to model C2 
and E2, the slope of the column density profile is nearly 
constant through entire region in model D2. 

FigurefTUlshows the magnetic field strength and plasma 
beta, which is defined as the ratio of thermal pressure 
to magnetic pressure, in these models. Although these 
three models have the similar magnetic field strengths of 
~ 150 fiG at cloud tip, the magnetic field strength de- 
creases more quickly with moving to tail side in the model 
with larger value of a. The distribution of plasma beta 
in the cloud qualitatively changes whether a exceeds 0.5 
or not. As can be seen from equation (0 and (JSJ, the 
thermal pressure and the magnetic pressure are propor- 
tional to p and p 2a , respectively. When a is larger than 
0.5, magnetic pressure increase faster than thermal pres- 
sure with increasing the density, and magnetic pressure 
become dominant at dense region. When a is smaller 
than 0.5, by contrast, thermal pressure become domi- 
nant at dense region. In model C2 in which a = 0.25, 
the plasma beta is larger than unity at dense head re- 
gion and decreases with moving to tail side. Thus, the 
thermal pressure is dominant at the head region; the 
magnetic pressure is dominant at tail side. In model 
E2 in which a = 0.75, opposite trend is observed. The 
magnetic pressure is dominant at the head region; the 
thermal pressure is dominant at tail side. In model D2 
in which a = 0.5, the plasma beta is nearly constant at 
entire cloud, because the magnetic pressure and the ther- 
mal pressure increase at the same rate with increasing of 
the density. 



4.3. Dependence on Initial Magnetic Field Strength 

In this subsection, we describe general trends which 
arise when initial magnetic field strength Bo is changed. 
In Figurc [TT1 and [T"2l maximum values of the number den- 
sity and the magnetic field strength are plotted as func- 
tion of initial magnetic field strength Bo, respectively. 
As shown in Figure [9] and 1 1 01 the number density and 
the magnetic field strength reach their maximum values 
near cloud tip, when effects of FUV radiation are not 
significant. We focus on low plasma regime and consider 
these maximum values as values at cloud tips. 

We can obtain asymptotic behavior of number density 
in low plasma beta regime by neglecting thermal pres- 
sure. From equation ([7]) and ([9]), the number density is 
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Fig. 9. — The density, the column density, and the temperature distributions in model A (left), B (center), and, E2 (right). 
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Fig. 10. — The magnetic field strength and the ratio of thermal pressure to the magnetic pressure in model C2 (left), D2 (center), and 
E2 (right). 



TABLE 2 
Summary of numerical results 



Model 


Rc [pc] 


T Up * [K] 


,max 


[cm''- 1 } 


Nn 2 ,max c [cm -2 ] 




A/ Ci c [M ] 


A 


7.6 x 




19.5 


4.48 x 


10° 


1.26 x lO^- 5 




61.3 


li 


9.4 x 


10~ 2 


28.9 


4.19 x 


10 5 


1.27 x 10 23 




54.8 


CI 


7.0 x 


HT 2 


30.3 


2.61 x 


10 s 


5.65 x 10 22 


119 


24.1 


C2 


4.1 x 


lO" 2 


29.5 


1.07 x 


10 s 


1.61 x 10 22 


172 


5.59 


C3 


3.2 x 


10" 2 


27.2 


5.35 x 


10 4 


6.36 x 10 21 


192 


2.03 


CI 


2.7 x 


10~ 2 


24.1 


2.27 x 


If) 4 


2.15 x 10 21 


208 


0.653 


Dl 


8.2 x 


10~ 2 


30.3 


1.81 x 


10 s 


5.33 x 10 22 


133 


27.7 


D2 


7.1 x 


10" 2 


27.5 


5.75 x 


10 4 


1.81 x 10 22 


161 


10.4 


D3 


6.7 x 


io~ 2 


24.2 


2.39 x 


10 4 


7.17 x 10 21 


173 


4.20 


Dl 


6.5 x 


lO" 2 


21.2 


7.83 x 


10 3 


2.27 x 10 21 


178 


1.35 


El 


9.0 x 


10~ 2 


30.1 


1.42 x 


10 s 


5.17 x 10 22 


139 


28.4 


E2 


8.9 x 


10" 2 


26.9 


4.87 x 


10 4 


2.05 x 10 22 


155 


13.0 


E3 


8.8 x 


10" 2 


23.8 


2.08 x 


10 4 


8.77 x 10 21 


164 


5.81 


El 


8.7 x 


10~ 2 


20.9 


6.35 x 


10 3 


2.65 x 10 21 


168 


1.79 



a Temperature at cloud tip 
b Maximum value of number density of H2 
c Maximum value of column density of H2 
d Maximum value of magnetic field strength 
c Total cloud mass 
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expressed as 



n Q (SttPu 



,l/2a D-l/a 
/ D ' 



(37) 



where we neglect the first term in the right-hand side of 
equation ([7]). because we assume that the magnetic pres- 
sure is dominant. Total pressure Ptot at cloud tip is given 
by substituting 9 = into equation (p~3|) . Therefore, us- 
ing the relation of Fjjv — Fi/q, we obtain the number 
density at cloud tip as 



rhi- 



no 



( IQirixmYtCiFi 



q 



l/2c 



B, 



-1/c 



(38) 



This relation shows that the density at cloud tip is pro- 
portional to B in low plasma beta regime. The maxi- 
mum number densities of our numerical models approach 
lines representing this analytic asymptotic value as Bq 
increases (see Figure lllj) . We adopted minimum value 
of R c for models with each a to calculate q in equation 
(|38|) . because curvature radius at cloud tip R c decreases 
with increasing of Bo. We used the value of R c = 0.027 
pc for models with a = 0.25, R c = 0.065 pc for mod- 
els with a = 0.50, and R c = 0.087 pc for models with 
a = 0.75. The estimate value of maximum number den- 
sity of SFO 22 is also plotted in this figure. Dividing max- 
imum value of observed column density 1.89 x 10 22 cm~ 2 
by cloud width of 0.26 pc gives rough estimate value of 
n bs — 2.36 x 10 4 cm -3 . From this figure, it is inferred 
that SFO 22 had initial magnetic field strength of several 
to ^90 /iG, if value of a is from 0.25 to 0.75. 

We can also obtain asymptotic behavior of magnetic 
field strength in low plasma beta regime. Substituting 
equation (|2"B1 and into equation (0) gives the mag- 
netic field strength at cloud tip as 



B 



tip 



( 16-7T (imnCjFi 
\ 1 



1/2 



(39) 



This relation shows that magnetic field strength at cloud 
tip is function of the curvature radius of the cloud at tip 
R c and incident ionizing photon flux Fi, but indepen- 
dent of the initial magnetic field strength Bo and initial 
number density no- As the initial magnetic field strength 
i?o increases, the magnetic field strength at cloud tip ap- 
proaches this asymptotic value. The equation (I3TR) gives 
asymptotic value of 212 fiG for models with a = 0.25, 
180 /iG for models with a = 0.50, and 170 ^G for models 
with a = 0.75. The maximum magnetic field strengths of 
our numerical models asymptotically approach this value 
at large B (see Figure [12]) • 

5. DISCUSSION 

5.1. Comparison between Observations and Numerical 

Models 

We developed numerical model for photo-evaporating 
cloud which is in quasi-stationary equilibrium state, as- 
suming that pressure gradient of the cloud balances the 
inertia force caused by acceleration due to the back reac- 
tion of photoevaporation flow. As shown in l4.lt observed 
density structure of BRC SFO 22 can not be explained 
by the reference model (model A) in which either effects 
of strong FUV radiation or effects of magnetic field were 
not included. The column density of the reference model 
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Fig. 11. — The maximum number densities as function of initial 
magnetic field strength. The circles (colored blue in the online ver- 
sion), the diamonds (colored red in the online version), the pluses 
(colored green in the online version), and the asterisk (colored light 
blue in the online version) represent model C1-C4, model D1-D4, 
model E1-E4, and model B, respectiv ely. The solid lines represent 
analytic estimates given by equation (1381) . 

(A color version of this figure is available in the online journal.) 
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Fig. 12. — The maximum magnetic field strengths as function 
of initial magnetic field strength. The circles (colored blue in the 
online version), the diamonds (colored red in the online version), 
and the pluses (colored green in the online version) represent model 
C1-C4, model D1-D4, and model E1-E4, respectively. The solid 
lines represent analytic asymptotic values given by equation (1391) . 
(A color version of this figure is available in the online journal.) 



is one order of magnitude higher than that of observed 
cloud. Observed lower column density implies that the 
pressure of observed cloud is higher than that of the refer- 
ence model. There are two possible explanations for this 
large difference of density structures between observed 
cloud and the reference model. The one explanation is 
that heating due to strong FUV radiation from exciting 
star warms the cloud, hence enhanced thermal pressure 
makes the cloud reach equilibrium state with lower den- 
sity than the reference model. Another explanation is 
that additional pressure due to magnetic field makes the 
cloud reach equilibrium state with lower density than the 
reference model. 

Comparisons between the reference model and numer- 
ical model including effects of strong FUV radiation 
(model B) show that strong FUV radiation has little in- 
fluence on structure of the cloud. Although strong FUV 
radiation slightly reduces the density and the column 
density of the cloud, the column density is much higher 
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than that of observed cloud. Small difference from the 
reference model suggests that heating due to FUV does 
not affect structure of the cloud at least in SFO 22. To 
affect on structure of the cloud, the photoelectric heating 
needs to be dominant heating process, i.e. r pe /r cr > 1. 
For rough estimation, let us approximate the efficiency of 
photoelectric heating e expressed by equation ([3"0"]) by a 
constant value of 4.87 x 10~ 2 . From equation ([29]). pi]), 
and the relation of tjjv = 2.5Ay, the ratio of photoelec- 
tric heating rate to cosmic ray heating rate is written 
as 

-^~43G exp(-2.5A;). (40) 

1 rr 



before compression 



Therefore, visual extinction of the cloud along the direc- 
tion to exciting star is need to be smaller than critical 
value 

Aw . MgM, (41) 

Applying this equation for SFO 22, it gives the critical 
value Ay cri as 3.3. The total mass traced by 13 CO emis- 
sion and width of SFO 22 give rough estimate of visual 
extinction of SFO 22 as 



A 



V,SF022 



12.0 Mr. 



/imH7r(0.13 pc) 2 ' 



(42) 



It gives the value of 15.4. This larger value of Av^sfoii 
compared to Ay.cri indicates that the heating due to 
FUV radiation affects only near cloud surface. This dis- 
cussion based on rough estimation is consistent with our 
numerical results. 

Comparisons between the reference model and numer- 
ical model including effects of magnetic field show that 
magnetic field strongly affects structure of the cloud. 
From results of model C2, D2, and E2, initial magnetic 
field strength of 5 - 45 is required to reproduce ob- 
served column density of ~ 2.0 x 10 22 cm -2 depending 
on the value of a. On the other hand, magnetic field 
strengt hs in molecular cloud s are measured by Zeeman 
effects. iCrutcher et al.l (|2010f ) statistically analyzed sam- 
ples of clouds with Zeeman observation in order to in- 
fer the distribution of the total magnetic field strength 
in the samples. According to their analysis, molecular 
clouds with the density of ~ 10 3 cm -3 have magnetic 
field strength of a few to several tens of /j,G. This co- 
incidence of magnetic field strengths between numeri- 
cal models and observations shows that observed column 
density of SFO 22 are naturally explained by effects of 
magnetic field. 

5.2. Magnetic Field Configurations in BRCs 

As shown in 14.21 structure of BRC strongly depends 
on the value of a in equation ((9|). The slope of column 
density profile at head region becomes steeper with de- 
creasing a. The slope of column density profile at tail 
region shows opposite dependence on a. Since this pa- 
rameter represents how much magnetic field is trapped 
in the gas during the compression, our numerical results 
implies that direction of magnetic field affects evolution 
of BRCs. When magnetic field is parallel to the direc- 
tion of UV radiation, the value of a is close to 0. On the 
other hand, when magnetic field is perpendicular to the 
direction of UV radiation, the value of a is close to 1. 




after compression 




Fig. 13. — Schematic figures of magnetic field configuration in 
SFO 22 before undergoing compression by radiation-driven implo- 
sion (top) and after reaching the phase of quasi-stationary equilib- 
rium (bottom). 



Since results of model E2 in which a was set to be 0.75 
shows best agreement with observations, magnetic field 
in SFO 22 is thought to be near perpendicular to UV 
radiation. However, deviation from observations can be 
seen at tail region. The slope of column density of model 
E2 is steeper than that of observed cloud at tail region 
(see Figure [8]) . This difference of column density pro- 
files implies possibility that the value of a a t tail r egion 
is smaller than head region. iMackev fc Lim] (pOll per- 
formed three dimensional MHD simulations and found 
that for weak and medium magnetic field strengths an 
initially perpendicular field is swept into alignment with 
tail during dynamical evolution. Their results may ex- 
plain the reason why the value of a is small at tail re- 
gion. Figure [13] shows schematic figure of magnetic field 
configuration in SFO 22 suggested by comparisons be- 
tween our numerical results and observations. The pos- 
sible scenario is as follows. Magnetic field in SFO 22 
was initially perpendicular to UV radiation from excit- 
ing star. Then, compression due to radiation-driven im- 
plosion make magnetic field be close to parallel to UV 
radiation at tail region. 

Very few attempts have been made at observationa l 
study on mag netic field of BRCs. iSridharan etal] (|1996T ) 
performed optical polarimetry observations toward CG 
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22 and reported that magnetic field is parallel to its 
tail. Their rough estimate gives magnetic field strength 
of ~ 30/iG. It is comparable with magnetic field strength 
obtained by our numerical models ( see FigurellOj) . Other 
optical polarimetry observations by lBhattl ()1999f ) showed 
that magnetic field in CG 30-31 is found to be nearly 
perpendicular to the cometary tails. Observational stud- 
ies on relations between density structures and magnetic 
field configurations in BRCs are required to reveal effects 
of magnetic field on evolutions of BRCs. 



3. Magnetic field strength and direction of magnetic 
field strongly affect structures of BRCs. Numerical 
model with initial magnetic field strength of 5 fiG 
shows the best agreement with the observations. 
When magnetic field is nearly parallel to UV radi- 
ation from exciting star, the cloud has steep column 
density profile at head region and flat column den- 
sity profile at tail region. When magnetic field is 
nearly perpendicular to UV radiation, the column 
density profile shows opposite trend. 



CONCLUSIONS 



Usi: 



sing the Nobeyama 45 m telescope, we observed BRC 
SFO 22 in the 12 CO (J = 1-0), 13 CO (J = 1-0), and 
C 18 (J = 1-0) lines. Observed column density profiles 
were compared with those of numerical models for photo- 
evaporating cloud in quasi-stationary equilibrium state 
in order to investigate how magnetic field and heating 
due to strong FUV radiation from exciting star affect 
structures of BRCs. We summarize our main conclusions 
as follows: 

1. From our radio observations, the column density 
profiles of SFO 22 along the line to its exciting star 
are nearly flat with the column density of ~ 10 22 

_2 

cm . 

2. Strong FUV radiation from exciting star has little 
influence on structure of SFO 22. Although en- 
hanced thermal pressure due to strong FUV radi- 
ation slightly reduces the density of the cloud, its 
effects arc not enough to reproduce the observed 
density structure of SFO 22. 



In this paper we only focus on quasi-stationary equi- 
librium phase of radiation-driven implosion model, and 
we will discuss implosion phase in a subsequent paper. 
We also plan further study using MHD simulation to es- 
tablish more realistic evolutionary model of BRCs. 
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